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We report on the analysis and prototype-characterization of a dual-electrode electro-optic mod- 
ulator that can generate both amplitude and phase modulations with a selectable relative phase, 
termed a universally tunable modulator (UTM). All modulation states can be reached by tuning 
only the electrical inputs, facilitating real-time tuning, and the device is shown to have good sup- 
pression and stability properties. A mathematical analysis is presented, including the development 
of a geometric phase representation for modulation. The experimental characterization of the device 
shows that relative suppressions of 38 dB, 39 dB and 30 dB for phase, single-sideband and carrier- 
suppressed modulations, respectively, can be obtained, as well as showing the device is well-behaved 
when scanning continuously through the parameter space of modulations. Uses for the device are 
discussed, including the tuning of lock points in optical locking schemes, single sideband applications, 
modulation fast-switching applications, and applications requiring combined modulations. 

PACS numbers: 



I. INTRODUCTION 

Although electro-optic amplitude modulators and 
phase modulators are commonplace in modern optics lab- 
oratories, there is no single, commercially available device 
that produces controllable amplitude and phase modu- 
lations with complete variability. When required, vari- 
ous optical configurations have been assembled, such as 
phase and amplitude modulators in series 0, and two 
phase modulators used in a Mach-Zchnder interferomet- 
ric setup (the Mach-Zehnder modulator) 0], which are 
both capable of amplitude and phase modulation. 

The device discussed in this paper, called a universally 
tunable modulator (UTM), is a modification of a com- 
mercial electro-optic amplitude modulator. An ampli- 
tude modulator consists of two phase-modulation crystals 
aligned at right angles, and connected with opposite po- 
larities to a single electrical input. A UTM (see Fig.^a)) 
is identical, but with two separate inputs, one connected 
to each crystal. A prototype was constructed from a New 
Focus broadband amplitude modulator (Model 4104) in 
this way. We demonstrate here that this extra degree 
of freedom enables the production of phase as well as 
amplitude modulation. 

We use the term " universally tunable" to highlight the 
fact that the device is capable of both amplitude mod- 
ulation (AM) and phase modulation (PM), with full se- 
lectability of the amplitude of each and of the relative 
phase between the two. Choosing a figure of merit with 
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which to assess the device is difficult in general, and de- 
pends on the modulation requirements of the application 
in question. We address this issue here by identifying 
two characteristics that might be expected of the modu- 
lations states produced by the device: variability and pu- 
rity. The variable nature of the device is brought out by 
experiments which require real-time tuning of the mod- 
ulator across its parameter space. Other experiments 
rely on precise selection of a small number of operating 
points, where the purity of the modulations produced has 
a direct influence on the success of the experiment. 

One anticipated use for the UTM device is in optical 
feedback control of the resonance condition of a Fabry- 
Perot cavity, or the fringe condition of a Michelson in- 
terferometer. The RF modulation techniques used to 
lock these devices are Pound-Drever-Hall locking Q and 
Schnupp modulation locking respectively, and are 
both based on properties of the devices that convert (in- 
jected) PM to AM, which is then demodulated to give 
an error signal readout. In both cases, the default lock- 
ing point is at a turning point in the transmitted (or 
reflected) optical field intensity; injection into the device 
of additional AM (along with the default PM) causes 
the device to lock with an offset relative to the turn- 
ing point. Moreover, the quadrature of AM required for 
Pound-Drever-Hall locking is orthogonal to that required 
for Schnupp modulation locking, introducing the possi- 
bility of independently tuning two lock points at once in a 
coupled system. This RF offset locking application is our 
case study example to demonstrate the variability prop- 
erty of a prototype UTM; the technique may be useful for 
future large-scale gravitational wave detector configura- 
tions [f| , with the feature of facilitating real-time tuning 
of detector frequency responses 
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The UTM device is capable of several "pure" states, 
which have various applications as well as providing a 
natural way to test the purity property of the proto- 
type. Single-sideband modulation, for example, requires 
an equal combination of AM and PM, in quadrature. The 
subject has arisen a number of times with optic fibre tech- 
nology where, for example, a single-sideband-modulated 
signal is immune to fibre dispersion penalties |2j, and 
the technique has also been suggested for subcarrier- 
multiplexing systems 0- Single-sideband modulation 
has been achieved by other means: optically filtering out 
one sideband |Toj : cascaded amplitude and phase modu- 
lators 0; Mach-Zehnder modulators 0; and more com- 
plex arrangements |Tlj . While these methods have merit, 
we submit that the UTM is far simpler in design, with 
fewer degrees of freedom available to drift, and is at least 
comparable regarding suppression capabilities. 



The UTM can produce PM or AM states, with purity 
limited by the accuracy of polarisation optics (includ- 
ing the birefringent effects of the device itself). While 
these states are obtainable using off-the-shelf amplitude 
or phase modulators, applications in coherent state quan- 
tum cryptography require fast-switching between AM 
and PM pjj, for which the UTM is ideal. Other ap- 
plications, including some quantum communication pro- 
tocols [l3j, would also benefit from having easy access 
to a tuned, stationary combination of AM and PM. In 
addition, we show that the UTM (or indeed an ampli- 
tude modulator) can produce a carrier suppression state, 
where the output consists only of two modulation side- 
bands. 



We have developed a geometric phase sphere descrip- 
tion of modulation states, the modulation sphere formu- 
lation, in analogy with the Poincare sphere description 
for polarisation states 0]. We demonstrate the use of 
this formalism by calculating the transfer function (elec- 
trical to optical) of the UTM device, both using optical- 
field-phasors and using modulation sphere parameters, 
and discuss the pros and cons of the two mathematical 
descriptions. We have found that the modulation sphere 
representation is particularly useful for gaining a visual 
insight into the dynamics of the UTM. 



Section [H] gives the mathematical description of the 
UTM, in terms of optical- field-phasors fSection lll B|) . and 
the modulation sphere formulation fSection III Cfl . The 
transfer functions from these two subsections are derived 
in Appendices IA II and IA 21 respectively. Section II I II de- 
scribes a characterisation experiment conducted on a pro- 
totype UTM, with the experimental layout described in 
Section fill Bl and the results presented in Section Till CI 
The main results of the study are summarised in Section 

El 



II. THEORETICAL MODEL 

A. Overview 

Before we proceed, a simple analogy may help to clar- 
ify the inner workings of a UTM device, and lend physical 
insight. A Mach-Zehnder modulator involves splitting a 
beam into two parts, separately phase-modulating each 
part, then recombining the two parts on a beamsplitter. 
The UTM is physically equivalent to this, where the two 
interferometer paths are collinear but different in polar- 
isation. Each UTM crystal modulates one and only one 
polarisation component. The phase between the elec- 
tric fields in the two polarisation states we identify as 
the interferometric recombination phase. The birefrin- 
gent waveplates used to create the initial polarisation 
take the place of the input beamsplitter, and a polar- 
ising beamsplitter facilitates the output recombination. 
The recombination phase is adjustable (by modifying the 
polarisation state), and one can control the amplitudes 
and phases of the single-beam phase-modulations gen- 
erated by each of the two crystals. Through manipu- 
lation of these input parameters, the user has complete 
control over the interference condition of the two carrier 
beams, and separately over the interference of the lower 
and upper sidebands generated by the modulating crys- 
tals. This amounts to having complete freedom to choose 
any modulation state by appropriate choice of the input 
parameters. 

As an example, consider two identically phase- 
modulated beams, where the optical phase of one beam is 
advanced by 90° , and the modulation phase of the other 
beam is advanced by 90°. The result, upon interfering 
the two beams, is that one resultant sideband is exactly 
cancelled out and the other is additively reinforced, pro- 
ducing a single sideband state. One can intuitively arrive 
at all of the modulation states discussed in this paper 
by reasoning along these lines. In the next section, we 
quantify such reasoning into a mathematical theory of 
the UTM. 



B. UTM Transfer Function 

As stated earlier, a UTM (see Fig. Qla)) consists of 
two modulating crystals in series, with separate voltage 
sources, and with modulating axes at right angles to each 
other. The UTM is used by sending a laser beam of ellip- 
tical polarisation through the system, and then through 
a linearly polarising filter angled at 45° to the crystals' 
modulating axes. 

The choice of polarisation state of the incident beam 
has a crucial effect on the transfer characteristics of the 
system. We will restrict ourselves to a subset of polar- 
isation states: those where there is equal power in the 
polarisation components aligned with each modulating 
crystal axis. Equivalently, we require that the polarisa- 
tion is elliptical with major and minor axes aligned at 45° 
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FIG. 1: (a) UTM schematic. Two modulating crystals are posi- 
tioned in series with modulating axes at 90° . A laser beam of ellip- 
tical polarisation passes through each crystal and then through a 
vertically polarising element, (b) The subset of input beam polari- 
sation states considered here is all states such that there is an equal 
optical field amplitude in both the left- and right-diagonal compo- 
nents, a equals the phase between these components: shown are 
(7 = (vertical), tt/4, it/2 (circular), 37r/4, and it (horizontal), (c) 
Phasor diagram representing the transfer function given by Eq. |3] 
If the input electrical signals are represented as rotating phasors (<5i 
and 82), then the corresponding phasors that represent the output 
phase and amplitude modulations (P and A) are proportional to 
their sum and difference respectively. 



to both crystal axes. For concreteness we will choose the 
ellipse axes to be horizontal and vertical, and the crystal 
axes to be left and right diagonal (as in Fig. ^a)). We 
define the angle a as being the phase difference between 
left-diagonal and right-diagonal electric field components 
of the light exiting the crystals (see Fig^b)). The iden- 
tity 



tan 2 ( : 



Ph 

Pv 



(1) 



follows, where Ph,v are the powers in the horizontal and 
vertical polarisation components of the input laser beam. 
We introduce this degree of freedom with the application 
of gravitational wave detector locking schemes in mind, 
where laser power is at a premium: choosing a such that 
the beam is almost vertically polarised allows us to retain 
the majority of the input carrier power, while retaining 
full access to all modulation states, albeit with reduced 
amplitude. 

Note that, in the ideal situation where both crystals 
have identical refractive indices and lengths, the phase 
angle a (and hence the polarisation state) is the same 
before and after the modulator. However, a real device 



almost certainly does not have identical crystals so that 
the relative phase between left- and right-diagonal com- 
ponents may well be changed by the crystals, in which 
case a should describe the polarisation state of the light 
exiting the modulator, just before it reaches the linearly 
polarising filter. 

The inputs to the two crystal electrodes are sinusoidal 
voltages of frequency aj m ; we describe these as com- 



plex phasors <$,• 



(j = 1,2, and 



-1 so 



that the single-crystal phase modulations are given by 
$t{Sje luJmt } = 5jCOs(uj m t + <j>j), where 3? and 3 respec- 
tively return the real and imaginary components of a 
complex variable. 

We can similarly write the output phase and ampli- 
tude modulations as complex phasors P — P e^ p and 
A = A e % ^ A where the optical field amplitude exiting 
the UTM can be shown to be: 



^OUT 



£1 



IN 



cos(cr/2) + m{Pe luJmt } + 9?{ie^ mt } 



= Bin [cos(er/2) + iPcos(cu m t + <j> P ) + Acos(uj m t + 4>a) ](2) 

where is the input optical field amplitude. In Eq. [21 
we have assumed small single-crystal modulation depths 
(5j <C 1), and we have factored away the net optical 
phase dependence. 

The relationship between input and output parameters 
for the UTM device can be shown to be: 



P = \cos 
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Si - 82 



(3) 



(A proof of this transfer function is outlined in Appendix 
EH) _ 

An important observation from Eq.|21is that the phases 
of the AM and PM track the phases of the sum and 
difference of the electrical signals, respectively. Thus, 
an electronic oscillator for use in demodulation schemes 
is readily available, and can be calibrated once for all 
operating points. This phase-tracking property holds for 
any choice of a; the effect of a change of polarisation 
ellipticity is to change the overall quantity of PM or AM 
available (as well as changing the overall beam power in 
Eq.0). 

Returning briefly to the topic of conserving as much 
carrier power as possible for gravitational wave detector 
applications, we derive a useful rule-of-thumb for select- 
ing an appropriate amount of power to tap off. Consider 
Eq. |3 when we set Si = iSi such that Si + S2 and Si — 5% 
have the same magnitude (we do this in order to focus on 
the coefficients of the complex Sj phasors). The ratio of 
modulation power between the AM and PM components 
is I Jl| 2 /|P| 2 = tan 2 (cr/2) and, comparing this with Eq.^ 
we have: 
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So, the fraction of power tapped off is directly propor- 
tional to the fraction of modulation power available to 
be expressed as AM (rather than PM). Typically then, 
if 10% of the output power is tapped off, and the user 
wishes to alternately produce a pure PM state and then 
a pure AM state (each of the same modulation power), 
then the AM state will require 9 times the input electrical 
power in compensation, compared to the PM state. 

It is apparent from Eqs. [5] and [3] that for a = it the 
carrier and the PM contribution vanish, leaving only the 
AM component. This is the carrier suppression operat- 
ing point, where only two sidebands remain. The mod- 
ulation ceases to be of the phase or amplitude type due 
to the absence of a carrier as a phase reference, though 
the modulation retains its beat phase through the second 
harmonic in optical power. 

A rearrangement of Eq. [21 serves to clarify the nature 
of single sideband modulation: 



E, 



OUT 
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IN 



cos(ct/2) + - (A + iP) 



iP 



(5) 



where the e U)m terms represent optical sidebands at 
±w TO relative to the optical carrier frequency, and an 
asterisk denotes a complex conjugate. It is clear then 
that choosing A = ±iP eliminates one or other of the 
sidebands. In other words, a single sideband state corre- 
sponds to having equal amounts of PM and AM where 
the two modulations are in quadrature. 



C. Modulation Sphere Formulation 

We now define a geometric phase representation of the 
modulation parameters of the system to further clarify 
the transfer characteristics of the UTM device. This is 
analogous to using Stokes parameters and the Poincare 
sphere representation for optical polarisation states 0] . 
In this case, we transform from the complex quantities 
P and A to a set of (real) coordinates (Mi,M 2 ,M 3 ) in 
M-space via the transformation: 

Mi = P 2 - A 2 

M 2 = 2PAcos(<j)p - 4> A ) = 2U{PA*} 

M 3 = 2PAsm{(j)p - <j) A ) = 23{Pl*} (6) 

A single point in M-space represents a distinct modula- 
tion state. In particular, the M-space representation sup- 
presses the common (beat) phase of the modulation state; 
the equations are defined in terms of the relative phase 
difference between AM and PM contributions. This is 
a useful simplification in that now every physically dis- 
tinct modulation corresponds to one point and one only. 
On the other hand, one cannot use the representation 
for demodulation-phase-tracking calculations, when the 
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FIG. 2: (a) Q-spacc diagram. Axis labels are visual reminders 
of relative amplitudes and phases of input electrical signals (try 
to visualise the display of a cathode ray oscilloscope run in X-Y 
mode with 8± and 82 as the X and Y inputs), (b) M-space diagram. 
Axis labels are optical amplitude phasors whose endpoints oscillate 
at the modulation frequency. For a = n/2, there is a graphical 
transfer correspondence between (a) and (b) according to Eq. [5] 
(c) In general, the transfer function of the UTM, Eq. [S] takes a 
sphere of constant input power in Q-space and returns an ellipse in 
M-space. A cross section is shown for a selection of values for a. 



renewed by the presence of an external (electrical) phase 
reference. 
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TABLE I: Details of significant operating points. The modulation parameters, expressed both in terms of phasors and M-space 
parameters, are given for each point. The electrical inputs required to produce these operating points are also given (both in 
terms of phasors and Q-space parameters), and these inputs vary with the polarisation parameter a. 



The three M-space parameters can be interpreted as 
follows: Mi measures the extent to which one kind of 
modulation (AM or PM) dominates over the other; M 2 
measures the degree to which the present PM and AM 
components are correlated or anti-correlated in phase; 
and M3 measures the degree to which the present PM and 
AM components are in quadrature phase to each other, 
and thus also measures the extent to which one frequency 
sideband has more power than the other. We also define 
the modulation power by Mo = P 2 + A 2 , from whence 
M 2 = M\ + Ml + Mf follows. Therefore, a surface of 
constant Mq is a sphere in M-space (see Fig. H^b)). 

We attach particular significance to modulation states 
whose M-space coordinates lie on one cardinal axis. The 
point (+Mo,0, 0) represents pure PM, and the point 
(-M , 0,0) represents pure AM. The points (0,±M ,0) 
represent correlated and anti-correlated PM and AM, 
and (0,0,±Aio) represent upper and lower single side- 
band states. Details of these states are summarised in 
Table I. These states, particularly those on the Mi and 
M3 axes, involve the suppression of a particular signal, 
and so provide a natural way to test the precision of a 
real device. 

Before we proceed, it is convenient to define an anal- 
ogous Q-space representation, (Qi,Q 2 ,Q3), for the elec- 
trical input parameters, via: 



Qi = 6 2 -S 2 

Q 2 = 25x8 2 cos,{(t>i-<t)2) = 2^t{8i5* 2 } 

Q 3 = 25 1 S 2 sm(4> 1 -h) = ^{5 1 S^} (7) 



where Qo = 8\ + S 2 — y/Q 2 + Q 2 +Q 2 is the sum input 
electrical power. Physical interpretations for these pa- 
rameters correspond with the M-space equivalents (see 
Fig. |2Ja) for a diagrammatic description). 

Now, we can rewrite the UTM transfer function, Eq.[31 



in terms of the M-space and Q-space parameters: 
M = J (Qo + Q 2 cos(ct)) 
Mi = i (Qocos(cr) + Q 2 ) 
M 2 = isin(cj)g 1 

M 3 = -hm(a)Q 3 (8) 
(A proof of this transfer function is outlined in Appendix 

E2) 

Since the modulation power Mo depends partly on 
Q2 (and hence depends on the detail of the input sig- 
nals, not just the overall input power Qo), a sphere 
of constant electrical input power in Q-space does not, 
in general, map to a sphere of constant modulation 
power in M-space. In fact, Eq. [8] describes a trans- 
fer function from a sphere in Q-space to an ellipsoid 
in M-space. It can be shown that the ellipsoid is cen- 
tered at (M 1 ,M 2 ,M 3 ) = (Q cos(cr)/4, 0, 0), has radii 
[Qo/4, Qosin(cr)/4, Q sin((j)/4], and has an eccentricity 
of e = |cos(ct)|. The ellipse is always prolate with long 
axis aligned with the Mi-axis, and always has one fo- 
cus at the M-space origin. The orthogonality of axes 
is preserved in the transformation, which consists only 
of a translation (the ellipsoid is not origin-centered) and 
a dilation (the ellipsoid is squashed in the M 2 and M 3 
directions). The centrepoint and proportions of the el- 
lipsoid are parametrised by the input light's polarisation 
parameter, a (see Fig. |2Jc)). 

A special case occurs for o~ = n/2 (circular input po- 
larisation), when Eq. |S] reduces to: 

M = \Qo 

M x = \q 2 

M 2 = ^Qi 

M 3 = -jQ 3 (9) 

Here, a sphere in Q-space maps to a sphere in M-space. 
This case is particularly instructive as Figs. and|2Ib) 
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take on a new significance: there is now a direct graphical 
correspondence between the two spheres, in accordance 
with Eq. So, as an example of using these spheres 
as a visual tool, we can see "at a glance" that the PM 
operating point is obtained by running both crystals with 
equal in-phase electrical inputs, or that a single sideband 
is obtained by running both crystals with equal inputs in 
quadrature. 



III. EXPERIMENTAL DEMONSTRATION 

A. Approach 

The following characterisation experiment was de- 
signed to measure the amplitude and phase response of 
a prototype UTM device with respect to both the AM 
and PM output states. In particular, we focus here on 
two kinds of measurement: those that measure the vari- 
ability of the UTM device, and those that measure the 
purity of the modulations that the UTM device can pro- 
duce. By variability, we refer to the capability of the 
UTM prototype to tune around the parameter space of 
modulations, or "dial up" a particular modulation, in 
a predictable, theory-matching fashion and with a rea- 
sonable level of precision. By purity, we speak of the 
UTM prototype's ability to accurately attain a particu- 
lar modulation (especially those that involve suppression 
of a frequency line) and to hold that modulation indefi- 
nitely without drifting. 



B. Experimental Layout 

The experiment is shown in Fig. [21 The polarisation of 
a laser source was prepared with a half- and quarter- wave 
plate in series, and the laser source was passed through 
the UTM, which we operated at w m /2-7r = 5 MHz. A 
polarising beam splitter completed the process; the ver- 
tically polarized output carried the modulation described 
in the theory above. 

The horizontally polarised output was used to keep 
track of the polarisation state, and hence to measure the 
value of a. The polarisation state used corresponded to 
a = 75° (near-circular polarisation with the vertical com- 
ponent stronger), for experimental convenience. 

Regarding the modulator itself, the original New Focus 
amplitude modulator (Model 4104) carried the specifica- 
tion: max V„ = 300 V @ 1.06 fjtm. While we did not 
explicitly measure for the prototype UTM (after the 
modification from the original AM device), we point out 
that it should still be of the same order of magnitude. 
Therefore, since the input voltages used did not exceed 
10 V, the following results are in the small modulation 
depth limit as described in Section HlH Bl 

A heterodyne detection scheme formed part of the 
measurement apparatus, where a shunted beam was 
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FIG. 3: UTM heterodyne characterisation experiment layout. A 
spectrum analyser provides frequency line data (of the direct AM 
beat, and heterodyned copies of sidebands) for establishing the pu- 
rity of the UTM , while two parallel double-demodulation circuits 
provide the complete set of modulation measurables for testing the 
variability of the UTM. BS: dielectric beam splitter; A/2 half-wave 
plate; A/4 quarter-wave plate; PBS: polarising beam splitter; PZT: 
piezo-electric actuator; AOM: acousto-optic modulator; FB: feed- 
back servo and HV amp; 90°: electronic phase-shifter; DC/HF: 
bias-T; SA: spectrum analyser. The 5 MHz signal generators are 
electronically phase-locked. 



frequency-shifted by cj^/27t = 80 MHz with an acousto- 
optic modulator (AOM). If this heterodyne oscillator 
beam is represented by £het = r yEi-^e luJht with 7 <§; 1 
and is interfered with the modulated beam from Eq. [2J 
then the detected power is given by: 



Pdet 


= Ef N [P DC + P Wm + P Uh + P UK -w m 






= cos (-); 






= 2cos(|)3?{ie^ mt }; 




p 


= 2cos(|)5R{ 7e ^ t }; 





Put 



p„ 



25? 
25R 



j(A + iP) 
7(i- ip) e l (^+^")i 



(10) 



where the power components are split up according to 
their respective frequencies. P u!m represents the measur- 
able beat due to the presence of AM, P [AJh is the beat 
between the carrier and the heterodyne oscillator, and 
Puih-Urn and P Uh + Um represent heterodyned copies of the 
two modulation sidebands. A spectrum analyser was 
used to monitor the strength of these frequency lines. 
The use of the spectrum analyser facilitated the assess- 
ment of the purity of states produced by the UTM pro- 
totype, by observing operating points that involved sup- 
pression of one of these frequency lines. 

An alternative set of measurables was obtained by elec- 
tronically mixing the heterodyne frequencies down to 
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baseband via a double- demodulation scheme. This gave a 
DC readout of the PM and AM amplitudes at a selected 
beat phase. A complete description of the modulation 
strengths and quadratures was obtained in this way, and 
was hence useful to characterise the variability of the 
UTM prototype. 

The method of signal extraction via double demodula- 
tion is best understood by reworking the last two com- 
ponents of Eg. HUIto: 



p 



P,„ 



2cos(-) p{Pe^™ t }S{7e l " ftt } 

+'0t{Ae luJmt }^{ 1 e luJht }\ll) 



This signal is mixed down to base band by demodulating 
at a;/j/27r = 80 MHz and then uj m /2n = 5 MHz in series. 
The output of a demodulation is sensitive to the relative 
phases of the signals being mixed: if they are exactly 
in-phase or out-of-phase, the output will be maximally 
positive or negative; if they are in quadrature, the out- 
put will be zero. Eq. ^2 shows that the u>h oscillator 
components of the AM and PM portions are orthogonal 
(one is the real component, the other is imaginary) so 
that the appropriate choice of electrical oscillator phase 
can force the readout of PM only, or AM only, or some 
linear combination of the two. Similarly, the demodu- 
lation quadrature of the uj m stage determines the beat 
phase that the output signal is projected onto. 

The output DC component of the first stage of demod- 
ulation (extracted with a bias-T component) was used 
as an error signal for locking the optical recombination 
phase of the heterodyne. This error signal varies sinu- 
soidally with respect to the optical recombination phase, 
with a zero crossing where the optical heterodyne beat 
{Pu> h ) and the electronic demodulation oscillator are in 
quadrature. In other words, the feedback loop "locks 
out" the heterodyne beat. Comparing Ea.llOland Eg. 1111 
we see that the AM term has the same phase as the P UJh 
term (that is, they both contain the terms 3t{je lulht }), 
so that the feedback loop also locks out the AM com- 
ponent of modulation, leaving only the PM component 
(which has a Q^e 1 "' 1 *} term). This is an important part 
of the process as, without a locking loop, the demodula- 
tion phase of the circuit would be uncharacterised, and 
would also be free to drift. 

In the experiment, two separate double demodulation 
schemes were used in order to measure PM and AM si- 
multaneously; the u>h signal was split with a 90° elec- 
tronic splitter to ensure that the two double demodu- 
lation circuits scanned orthogonal modulations. Hence, 
when the heterodyne locking loop was connected using 
a signal from one double-demodulation circuit, this en- 
sured that that particular circuit was sensitive to PM 
and that the other was sensitive to AM. 

Regarding the second stage of demodulation, electron- 
ically phase locking the signal generators was sufficient to 
time-stabilise the demodulation quadrature, and a known 
modulation was used for calibration. 
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FIG. 4: Spectrum analyser traces of heterodyne measurements 
of four "pure" operating points. Frequencies of interest are 5 MHz 
(direct AM beat), 80 MHz (heterodyne-carrier beat), and 75 MHz 
and 85 MHz (heterodyne-sideband beats). A solid line indicates 
data averaged over a few seconds, a dotted line indicates max-hold 
data acquired over approximately an hour, and arrows indicate 
lines that have been suppressed. The plots are: (a) upper single- 
sideband (lower sideband suppressed by 34.9 dB); (b) lower single- 
sideband (upper sideband suppressed by 39.5 dB); (c) pure PM 
(AM beat suppressed by around 38 dB); (d) carrier suppression 
(heterodyne beat suppressed by around 30 dB). 



C. Results 

Fig. 0] shows spectrum analyser traces for four oper- 
ating points, to demonstrate the suppression capabilities 
of the device. The suppression factors are 34.9 dB and 
39.5 dB for the left- and right-hand sidebands (Fig. 0Ja) 
and Fig.^fb)). The figures shown suppress the frequency 
lines down to the electronic noise floor, so that higher 
suppression factors may be possible. A more detailed 
trace exhibiting sideband suppression is given in Fig. 
with 35.2 dB relative suppression recorded. 

We measured around 38 dB suppression of the AM 
beat compared with a similar-input-power pure- AM state 
(Fig. 0Jc)). In this and the other diagrams in Fig. 01 we 
have shown " max-hold" data, demonstrating that the de- 
vice is highly stable, maintaining these operating points 
without significant drift on the timcscalc of hours. 

In Fig. 01d), the carrier is suppressed by selecting the 
AM operating point and setting a = n (horizontal po- 
larisation). The heterodyne measurement of the carrier 
is down by around 30 dB, from which we can infer that 
the carrier power (which goes as the square of the hetero- 
dyne measurement) is down by 60 dB. This heterodyne 
measurement is further supported by observing that the 
first and second harmonics of direct AM beat are approx- 
imately equal, which is consistent with the fact that the 
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FIG. 5: A more detailed spectrum analyser trace of sideband 
suppression, with 35.2 dB difference between sidebands. Frequency 
lines at 72.1 MHz, 86.1 MHz and 87.7 MHz are from radio inter- 
ference with the electronic equipment. 



carrier (as measured by the heterodyne) is approximately 
6 dB weaker than the sidebands. 

We were especially careful to validate this heterodyne 
measurement of the carrier, because the directly detected 
power only dropped by 40 dB. The reason for this is that 
the majority of the residual carrier power was now in a 
higher-order, odd spatial mode (easily verified by looking 
at the intensity profile of the beam), which did not inter- 
fere efficiently with the heterodyne oscillator beam. The 
higher-order modes were thought to result from a spa- 
tially non-uniform polarisation of the light exiting the 
UTM which, when subsequently passed through a polar- 
ising beam splitter, produced modes reflecting the sym- 
metry of the modulator. As such, for applications where 
spatial mode interference is important, the employment 
of a mode-cleaner cavity with free spectral range equal to 
the modulation frequency (or one of its integer divisors) 
should solve the problem. 

In an attempt to confirm the variability of the UTM 
prototype, we endeavoured to map significant paths in 
the modulation parameter space. Fig. Eland Fig. [7| show 
the result of sweeping through M-space in two cardi- 
nal directions: varying the relative phase between the 
two (constant amplitude) electrical signals; and varying 
the amplitude of one electrical signal while keeping the 
other signal amplitude, and the relative phase, constant. 
For matters of calibration, the UTM was initially set 
to the PM operating point (this initial state was used 
as the reference in labelling "in-phase" and "quadra- 
ture" components), and the relative gains and phases 
of the two double demodulation circuits were measured 
and factored out via correspondence with spectrum anal- 
yser data. The input polarisation was set to be elliptical 
with a vertical component slightly larger than the hori- 
zontal component. The results show that changing the 
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FIG. 6: (a) Modulation ellipse showing trajectory of measurement 
sweep through parameter space, where the phase difference between 
electrical signals, <f> = <f>i — <f>2 , is varied through 360° . (The ellipse 
has an (exaggerated) a of 53°.) (b) Double-demodulation mea- 
surements of PM (squares and diamonds for in-phase and quadra- 
ture components respectively) and AM (pluses and crosses for in- 
phase and quadrature components respectively), with correspond- 
ing theoretical predictions (solid lines). The majority of modula- 
tion present is either PM or quadrature AM as predicted. This 
plot corresponds to a 75°; a value derived by measuring, from 
corresponding spectrum analyser data, the phase <j> at which the 
single sideband operating points occur. Using a result from Table 
I, the phase (f> between the pure AM operating point and either of 
the single sideband operating points, is precisely equal to a. 



phase between the two electrical signals produces AM 
in quadrature to the original PM, whereas changing the 
amplitude of one electrical signal generates AM in-phase 
with the original PM. The data points are generally in 
good agreement with the theoretical predictions. Some 
small systematic errors are apparent, and are thought 
to be associated with unmatched electrical and optical 
impedances between the two crystals, and possibly to 
do with slightly unmatched optical power levels probing 
each crystal (ie a violation of our original assumption re- 
garding polarisation states) . Overall, the device is shown 
to be highly predictable, and that one can "dial-up" a 
particular modulation state on call, having calibrated the 
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FIG. 7: (a) Modulation ellipse showing trajectory of measurement 
sweep through parameter space, where the strength of one electri- 
cal signal, 82, is varied from +5.1 V to —5.1 V. The sweep involves 
varying the overall input electrical power, so the modulation tra- 
jectory does not stay on the ellipse surface, and instead traces out 
a parabolic curve. (The ellipse has an (exaggerated) a of « 53°.) 
(b) Double-demodulation measurements of PM (squares and dia- 
monds for in-phase and quadrature components respectively) and 
AM (pluses and crosses for in-phase and quadrature components re- 
spectively), with corresponding theoretical predictions (solid lines). 
There is diverging agreement in the PM which may be caused by 
imperfectly matched electrical impedances for the two crystals. As 
in Fig. |n] the location of the correlated PM & AM point is shifted 
toward the pure AM point by about 0.7 V out of 5.1 V and, refer- 
ring to Table ITTT?1 this corresponds to a value of a ss 75°. 
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FIG. 8: (a) Modulation ellipse showing short-range trajectories 
designed to map local region near PM operating point. The elec- 
trical parameters <f> = <f>i — <p2 and 82 were varied to achieve these 
results, with units worth 5° and 0.44 V respectively, (b), (c), (d) 
and (e) show double-demodulation measurements of PM (squares 
and diamonds for in-phase and quadrature components respec- 
tively) and AM (pluses and crosses for in-phase and quadrature 
components respectively), with corresponding theoretical predic- 
tions (solid lines). Parameters varied were (b) 0, (c) 82, (d) <f> and 
82 together with the same polarity, and (e) </> and 82 together with 
opposite polarities. Note: the in-phase PM data points (squares) 
have been scaled down by a factor of 10 to fit in the diagram. Also, 
the systematic error in the quadrature PM data is probably due 
to an electronic phase drift between signal generators, causing a 
small amount of in-phase PM data (which is an order of magni- 
tude stronger) to couple across. 



device initially. 

It is instructive to notice that the cardinal operating 
points are not evenly spaced when sweeping across the 
parameter space. In particular, we see that the two 
single-sideband operating points are shifted toward the 
AM operating point in Fig. El and the correlated PM & 
AM operating point is similarly shifted in Fig. [7| rela- 
tive to the perpendicular. This is a product of the non- 
circular polarisation of light input into the modulator de- 
vice, and can be understood by reviewing the modulation 
ellipses in the two figures; the Mi and M3 axes intercept 
the ellipse at points that are closer to the AM operating 



point than the PM. In fact, this provides us with a means 
of calibration of the relative transfer- function-amplitudes 
of the PM and AM double-demodulation circuits: they 
are set by assuring that (in Fig.[5J) the PM and AM com- 
ponents are equal at the phase <j) where the single side- 
band is known to occur (information which is obtained 
via comparison with spectrum analyser data). Also, by 
using the identities in Table I, we can derive a value of 
the polarisation parameter: a rs 75°. 

Fig. [8] shows scans of the local region near the PM op- 
erating point, with deviations in four directions. These 
data were taken with an offset locking application in 
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mind, as discussed earlier. The in-phase PM component 
(in-phase by definition) is by far the dominant signal, 
and has been scaled down by a factor of 10 to fit in the 
graphs. As expected, the UTM device can produce AM 
that is in-phase with the present PM, or AM that is in 
quadrature with the present PM, by changing the relative 
electrical signal amplitudes or phases respectively. In ad- 
dition, both of these parameters can be varied together 
or oppositely, to give independent control over the two 
quadratures of AM resulting. We note that in the case 
of producing both in-phase and quadrature AM, a small 
component of quadrature PM appears (in other words, 
the phase of the PM changes, relative to that of the pure 
PM point used for the initial calibration). In actual fact, 
the quadrature PM data has a noticeably larger system- 
atic error than the other three signals. This is most likely 
due to pollution from its in-phase counterpart, caused by 
the demodulation oscillator's phase drifting marginally 
(this was observed to happen even in spite of the elec- 
tronic phase locking between signal generators). 

A number of experimental difficulties deserve mention- 
ing. As described in the context of the carrier suppres- 
sion results, the UTM produced a spatially varying po- 
larisation state, which produced a small percentage of 
higher-order spatial modes upon selecting out a polar- 
isation component (of order 1% of the overall power). 
This interfered with our ability to directly measure the 
polarisation, which we did by measuring the overall de- 
tected power while rotating a diagnostic half-wave plate 
placed before the PBS. In this way, we measured a value 
of a » 70° for the results shown in Figs. [HI to [SJ as 
compared to a ~ 75° measured by inference from spec- 
trum analyser data. We assign a reasonably large error 
to the direct measurement value due to the presence of 
the higher-order spatial modes, and we favour measuring 
cr by inference from the data. 

The overall optical path lengths of the modulat- 
ing crystals were found to change significantly as they 
warmed up after starting the laser. The dual crystal de- 
sign of the device goes a long way toward minimising 
this problem, and we have seen that the device is stable 
once it has warmed up. However, long term drift of the 
crystal lengths is possible, having a direct effect on the 
polarisation state leaving the device. We propose that, 
in circumstances where this becomes a problem, a feed- 
back loop be employed to lock the polarisation state. A 
possible scheme would see the power level out of one PBS 
port monitored and used as a feedback signal (minus an 
offset equal to the desired power level) to the DC optical 
path length of one of the crystals, thus compensating for 
a mechanical path length change by feeding back to the 
refractive index. 

Finally, the electrical impedances of the two crystals 
were not well matched for two reasons: the electron- 
ics were not identical; and the crystals themselves had 
"good" and "bad" spots which generated varying levels 
of modulation. In general terms, careful alignment can 
largely overcome this problem and return the two crystals 



to an equal footing. The differing electrical impedances 
complicate the issue of generating a local oscillator with 
the phase-tracking property discussed earlier, since the 
voltages Vi,2 and modulation depths are not then 
related by the same factor. We submit that the device 
tested here is merely a prototype, and that careful man- 
agement will suffice to deal with these issues as the need 
arises. 



IV. CONCLUSION 

We presented a thorough investigation, both theoret- 
ical and experimental, of a prototype Universally Tun- 
able Modulator (UTM). The electrical-to-optical trans- 
fer function of the device was derived, both in terms of 
electrical and optical phasor notation, and using a geo- 
metric phase "Modulation Sphere" representation. Both 
pictures were shown to have merit, and a set of cardi- 
nal modulations were described in each notation. We 
reported on an experiment to characterise the prototype 
UTM, which involved using dual double-demodulation 
circuits to measure both AM and PM components si- 
multaneously. Data sets were obtained and analysed to 
illustrate the variability and purity characteristics of the 
device; the device was shown to be highly predictable 
and capable of highly pure states. Applications for the 
UTM were discussed. 
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APPENDIX A: APPENDICES 
1. Derivation of UTM transfer function 

This is an outline of the derivation of Eq. [3J First, 
we find the (more general) transfer function of the UTM 
when we allow the input beam to have any polarisation. 
At the end, we will simplify to the subcase described in 
the main text. 

The input beam's polarisation will be characterised by 
two electric field phasor components, L = Le l<TL and 
R = Re ltTR , as defined in a set of left-diagonal and right- 
diagonal spatial coordinate axes (with unit vectors L and 
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R respectively). The electric field exiting the UTM can 
be written as a vector (to include polarisation informa- 
tion) as: 



Ei 



ExitingUTM 



(Al) 

(From hereon in, we will suppress the e lu)t term for 
the sake of brevity.) Upon passing through the verti- 
cally aligned linear polariser (equivalent to taking a dot 
product with the vector [L + R)/y/2), and expanding 
the complex exponentials to first order (hence assuming 
|<Si,2 1 "C 1), the vertical electric field amplitude becomes: 



&r — <Jl as the phase between the two polarisation com- 
ponents, then the equations reduce to: 

£out = ^incos (I) + &{Ae* u '" t } + miPe^} (A7) 



and 



- E m ( a 
P = —cos 



A 



(§)(i + 4) 

Si - S 2 ) • (A8) 



Eqs. 13 and [3] in the text are the same as these, but 
with a dimensionless definition of P and A. 



(A2) 

Next, we collect DC terms and oscillating terms, and 
factor out the overall optical phase: 



Eqvt — 



L + R 

~7T 



i(RL* + R 2 ) 
\L + R\ 2 



\L + R\ 2 



n{6 2 e luJmt } 



(A3) 



Now, wc work to separate real and imaginary compo- 
nents for the oscillating terms: 



K 



OUT = 



L + R 



1 + S W^e^' 

V \L + R\ 2 



( nm& + s 2 )+L% + R%) \i 

v l^+^l 2 7 J 

These real and imaginary oscillating terms correspond 
to AM and PM respectively, so we parameterize via: 

£ OUT = |L ^ 1 + U{Ae iUmt } + iSt{Pe iUmt } (A5) 



with 



P = 3? 



V2 



/ 3?{i?L*}(5i + S 2 ) + EH X + R 2 5 2 ) ' 



.4 



3J 



( %{RL*}(6 1 -~5 2 y 
{ V2\L + R\ 



R\ 



(A6) 



where the net optical phase shift has been discarded. Eqs. 
A5l and lA"6l constitute the general UTM transfer function 
for arbitrary input polarisation. Note that the PM com- 
ponent does not have the property of tracking the phase 
of the sum of the input electrical signals, since it gen- 
erally depends on these inputs in different proportions. 
This is one of the reasons why we chose to restrict the 
polarisation to a subset. 

If we choose said subset, L = R = Sin/v^ (equal 
power in the two polarisation axes), and define a = 



2. Derivation of UTM transfer function for 
Q-space and M-space parameters 

This is an outline of the derivation of Eq. [S] As with 
Appendix I A II we derive the Q-space to M-space trans- 
fer function of the UTM for any input light polarisa- 
tion state. At the end, we reduce the equations to the 
case where the allowed polarisations are restricted. The 
derivation consists of transforming parameters P and A 
to M-space parameters (Mi, M 2 , M3) using Eq. El and 
transforming parameters <5i and S 2 to Q-space parameters 
(<3i,Q2,Q3) using Eq. [3 and hence converting Eg. IA6I 
from one set of coordinates to another. 

Firstly, it is convenient to write down a few interme- 
diate terms using Eq. IA6I as a starting point. We also 
convert to Q-space parameters (via Eq. [7J as we go. 



p2 = - = $t{RL*} 2 [Q Q + Q 2 ] + 2L 2 R 2 Q 2 



2\L + R\ 2 



i[(L 4 + i? 4 )Q + (£ 4 -i? 4 )Qi] 



2\L + R\ 2 

^{RL*}[(L 2 + i? 2 )Q + (L 2 - R 2 )Qi + (L 2 + R 2 )Q 2 ] 



2\L + R\< 



A9) 



,2 1. 7 Z{RL*} 2 [Qo-Q 2 ] 



A 2 = A* A = 



2\L + R\ 2 



(A10) 



ft{PA*} = 



ft{RL*}Z{RL*}Qi 



2\L + R\ 2 

&{RL*} [(L 2 - R 2 )Q + (L 2 + R^Q, - (L 2 - ffig^ ■ 
2|L + i?| 2 



%{PA*} = 



{RL*} n{RL*} + \{L 2 + R 2 ) 



2\L + R\< 



^A12) 



Things can be simplified somewhat by using a Stokes' 
parameter representation for the polarisation from this 
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point. Here, we (unusually) define the Stokes' parameters 
in terms of left-diagonal and right-diagonal electric field 
components (equivalent to rotating the usual X-Y axes 
anti-clockwise by 45°): 

5 = R 2 + L 2 

51 = R 2 -L 2 

5 2 = 2RLcos(a) = 2$t{RL*} 

5 3 = 2i? J Lsin(cr) = 2%{RL*} 

when the previous four equations become: 
1 



(A13) 
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(A14) 



These intermediate quantities are then converted to M- 
space parameters (via Eq. |BJ . We write the result here 
in matrix form: 
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(A15) 



Eq. IA15I is the general Q-space to M-space transfer 
function for the UTM with any input polarisation. It 
gives further insight into the manner in which choosing 
L 7^ R affects the properties of the system. The pa- 
rameters Mi and M 2 depend on both Q\ and Q 2 , which 
means that a sphere in Q-space will be distorted in the 
Afi-M 2 plane by the transformation. The resulting sur- 
face is an ellipsoid whose major axis is no longer collinear 
with the Mi-axis, but is at an angle subtended in the 
Mi-M 2 plane. This distortion certainly interferes with 
most of the favourable properties of the system. For ex- 
ample, none of the poles of the ellipse coincide with a 
coordinate axis, making any of the six operating points 
discussed more complicated to find. 



The restriction on polarisation states discussed in the 
text is equivalent to setting Si to zero (which in turn 
forces So = £f N , S 2 = £' I 2 N cos(cr) and S3 = £' I 2 N sin(cr)), 
when Ea. lA15l reduces to: 
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Once again, Eq. [5] from the text is the dimensionless 
equivalent of this, where -Ein is external to the definition 
of the M-space parameters. 
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